Ernst equation and spheroidal coordinates with a 
cosmological constant term 

o ' 

C.Charmousis 



(N 



O . LPT, Universite de Paris-Sud, 

Bat. 210, 91405 Orsay CEDEX, France 
LPT-0672 



o 



o 



E-mail: christos . charmousis@th . u-psud . f r 



cr 

Abstract. We discuss solution generating techniques treating stationary and axially 
symmetric metrics in the presence of a cosmological constant. Using the recently found 
extended form of Ernst's complex equation, which takes into account the cosmological constant 
term, we propose an extension of spheroidal coordinates adapted to asymptotically de-Sitter 
and anti de-Sitter static spacetimes. In the absence of a cosmological constant we show in 
addition that any higher dimensional metric parametrised by a single angular momentum can 
be given by a 4 dimensional solution and Weyl potentials parametrising the extra Killing 
directions. We explicitly show how a stationary, and a static axially symmetric spacetime 
solution in 4 dimensions, can be added together to give a 5 dimensional stationary and 



axisymmetric solution. 



^ , 1. Introduction 

In recent years there has been an increasing effort in finding exact solutions of higher 
dimensional gravity [1], [2], [3], [4]. In particular, since the pioneering work of Maldacena 
[5] bringing into perspective the adS/CFT correspondence some effort has been devoted to 
understanding the effect of the cosmological constant term in Einstein gravity (see for example 
[6] , [7] , [8] ) . Such studies have also been motivated recently from braneworld gravity. Using the 
adS/CFT correspondence in the context of braneworlds [9], intriguing relations between bulk 
higher dimensional black holes and their 4-dimensional quantum versions (see also [10]) have 
been put forward. If this correspondance were true, exact higher dimensional solutions solutions 
would be providing intriguing information about the quantum description of 4 dimensional 
black holes from higher dimensional classical solutions. However (and not-surprisingly) such 
solutions have been proven very difficult to find analytically, and in particular in the presence 
of a cosmological constant term, which is vital for an adS/CFT description. It seems that 
better understanding of solution generating techniques as well as the investigation of convenient 
coordinate systems, involving the cosmological constant term, are needed in order to tackle such 
problems were, even in the case of 4 dimensional general relativity, very little is known. 



In a recent paper [8], rotating spacetimes of axial symmetry were studied in the presence of a 
cosmological constant. Classical techniques, as that of Lewis-Papapetrou [15], were developed 
there to include the cosmological constant term. In this letter, we will firstly focus on the Ernst 
equation [16], as well as extend on a solution generating method developed in [8]. The Ernst 
equation, extended for a cosmological constant term [8], will permit us here to propose an 
extension of spheroidal coordinates to adS/dS static black holes. These coordinates, for A = 0, 
have been shown to be very useful in the study of stationary metrics. In particular Ernst [16] 
showed how one could generate rather simply Kerr's solution starting from Schwarzschild. The 
extension of asymptotically flat coordinate systems to asymptotically A / coordinates maybe 
very important in order to find novel stationary solutions such as the adS version of the black 
ring solution [2]. The latter solution generating method on the other hand will allow us, using 
Weyl's classical GR formalism, to extend 4 dimensional solutions to higher dimensional ones. 
In particular we will show how by literally adding together the potentials of 4 dimensional 
static and stationary metrics we can construct 5 dimensional stationary ones. 

In 4 dimensional general relativity Einstein's equations in the vacuum Rab = 0, guarantee 
that any locally static and axially-symmetric metric can be written as 

where the metric components depend on r and z. Since the field equation for a reads, 

Aa = (2) 

by a suitable 2 dimensional conformal coordinate transformation we can set a = R thus 
obtaining the Weyl form [12] 

ds 2 = -e 2X dt 2 + e" 2A [R 2 d(p 2 + e 2 *(dR 2 + dZ 2 )] , (3) 

where A, the Weyl potential, and x depend on R, Z and satisfy the field equations, 

di + ±d R + d 2 \\ = 0. (4) 

dnx = R [(<9rA) 2 - (d z A) 2 ] , dzx = 2Rd K Xd z X (5) 

Equation (4) is the linear Laplace equation written in three dimensional cylindrical coordinates. 
Given that (4) is a linear equation, a general solution can be found by simple separation of 
variables and imposing adequate asymptotic boundary conditions [13]. In a more pictorial 
manner, one can view Weyl potentials A, as Newtonian sources in three dimensions and can, 
according to (4) even superpose them creating new solutions from known ones. Once the 
Weyl potential A has been specified, one evaluates the x field by direct integration from (5). 
Equations (5) actually carry the full non-linearity of Einstein's equations. Weyl potentials are 
not unique and are associated to the patch of coordinates we are using. Flat space for example 
has Weyl potentials (modulo a constant) given by, 



A = 0, A = In R, A = \ ln^R 2 + Z 2 + Z) (6) 

where in particular the last one is adapted to an accelerating Rindler patch. A Schwarzschild 
black hole of mass M has Weyl potential given by 

x _ 1 , / R + +R,-2M \ m 



where R± = R 2 + (Z ± 2M) 2 . As we mentioned one can superpose black hole Weyl potentials 
obtaining multiple black hole solutions [14]. Typically when sources are superposed, conical 
singularities for R — ► appear, and are interpreted as struts holding, or strings pulling the 
sources apart thus keeping the attractive sources in a static equilibrium. 

In essence, Weyl components adapt the problem of finding solutions to a three dimensional 
flat coordinate system. This system is convenient for analysing the solutions in parallel with 
their Newtonian sources but is not always tailored to the solutions themselves. Quite often it 
is more suitable to adapt the coordinate system to the Weyl potential A rather than a. This is 
one of the ideas behind spheroidals although they were not initially introduced or defined this 
way. These coordinates were first discussed in the context of axial symmetric spacetimes by 
Zipoy 1 [19] (see also [3] for D > 4). Consider polar dike coordinates (u,ip) but with hyperbolae 
as radial functions, that is 

Z = cosh u cos tp , 

R = sinh u sin ip , (8) 

so that in the (R, Z) plane ip = const curves are hyperboloids and u = const are ellipsoids. On 
setting x = coshu and y = cos ip, the coordinate system becomes anew symmetric in x and y. 
Consider a Schwarzschild black hole: the standard metric, 



ds 2 = - 1 )dt 2 + ^ + r 2 dnjj (9) 



2 , „ 2M\ , 2 dr 2 2 ,^2 

- * —) '" ■ r ~JE 

can be rewritten in Weyl coordinates (r,z) of (1) setting r 2 /2M = cosh 2 (r/2) and 6 = z. The 
conformal transformation to (3) and (8) gives 

e 2A = (10) 
x + l y ' 

so that in spheroidals the Weyl potential for (9) is rather simple. We will come back to this 
point in a moment but first let us switch-on rotation. 

Lewis and Papapetrou [15] generalised the approach of Weyl to stationary and axisymmetric 
solutions in vacuum. After a conformal transformation, the metric takes the Lewis-Papapetrou 
form 

ds 2 = -e 2X (dt + Adtp) 2 + e' 2X [R 2 d(p 2 + e 2x (dR 2 + dZ 2 )] , (11) 

which differs from the static form by the additional component A = A(R,Z). Note that dt is 
no longer a static but rather a stationary locally timelike Killing vector field. For the metric 
(11), Ernst [16] pointed out an interesting reformulation of Einstein's equations for A and A, 
which read respectively 

9r (^<^) + dz {^dzA^j = 0, (a 2 + ±d R + a 2 ) A = |^ [(dnA) 2 + (d z A) 2 ] . 
, • ' ^ " - * (12) 

Indeed introduce an auxiliary field, u, defined by 

4A 

(-dzu, dRu;) = — (c\A, d z A) , (13) 



1 According to [19] such coordinates were used to describe the exact Newtonian gravitational field of the Earth. 



and the complex function 



then satisfies the differential equation 



£ = e 2X + iuj 



(14) 



(15) 



where ^ = (Or, <9z). This complex partial differential equation is known as the Ernst equation 
[16]. Its real and imaginary part are exactly (12). In this language, the Weyl potential A 
is simply given by the real part of the Ernst potential £ , whereas rotation is embodied by a 
non-trivial to (or A). 

Using the symmetries of complex functions, several methods have been proposed to obtain 
solutions of the Ernst equation (15) and hence to generate new solutions (see [11], [16], [17] 
and references within). An elegant example appeared in Ernst's original paper [16], namely a 
simple method to obtain the Kerr solution from the Schwarzschild solution using spheroidal 
coordinates. Indeed, consider the Mobius transform, 



£ = 



£ + 1" 



(16) 



defining a new potential £ which now solves the Ernst equation, 



2£* (V^ 



lei 5 



i 



(17) 



where a star denotes complex conjugation. 2 It follows immediately from (10) and (14) that 
£ = x for (9). Hence we notice that our transformed Ernst potential £ is now the new 'radial' 
coordinate x. This is in contrast to (1) where a = R and the Ernst potential is given by (7). 
In other words we have adapted the coordinate system to the real part of the black hole Ernst 
potential. We will be using this as our starting definition for extending spheroidal coordinates 
when we will switch on the cosmological constant term. Just in order to complete the discussion 
on the construction of Kerr, given the x <-> y symmetry, £ = y is also solution of (15), as is 
£ = xsini9 + iy costf. It turns out that this is nothing other than the Ernst potential of the 
Kerr black hole, where sini? = a/M is the ratio between the angular momentum parameter 
and the mass of the black hole [16]. 

When we consider a D-dimensional spacetime with a cosmological constant, a rotating 
metric of axial symmetry and with a single component angular momentum can be conveniently 
written in the form [8], 



ds 2 



2u - 

e a 



D-2(dr 2 + dz 2 ) + a°- 



' V2CE 



(D-2) 



*0 



-e~2(dt + Adip) 2 + e 



-SI 

2 d^ 



+ 



+ e 



2)(D-4) 



*0 



D-4 



(18) 



We have (D — 2) Killing vectors but as for D 
dt and &&. The fields a, v 



4, A = only two of them are not orthogonal, 



O and ^ again depend on r and z. These metric components 



2 In this representation of the potential, (17) is invariant under the complex transformation £ — > ^e 1 ^ for any 
phase ■& G R and stationary solutions can be easily generated from static ones (see [16] and [8] for more details). 



extend the Lewis-Papapetrou form of the previous section [8]. Indeed the field equations take 
the form 



Aa = -2ka— e 2v (19) 

= ^ • (ePa^A) (20) 

^■(ctfri) = 2ee n (^Ay (21) 

\?-(aVv^ = 0, n = 0...d-3 (22) 

= \ (*L + + ( A J 2 + E *?«. (« ~ «) ( 23 ) 



2^,„ — 



a a 



(24) 



If A 7^ 0, a is no longer a harmonic function. Furthermore, a fixes f from (19) which is no 
longer directly given by (23) as for A = 0. Pictorially a is a messenger component relating 
all the equations together. This is in contrast to the case when A = and equations (22) are 
decoupled from (20) and (21) (see the geometric interpretation in [8]). The components ^ 
and fi play a similar role to the Weyl potential A of (3). When D = 4 we have by definition 
^ = 0. It will be useful here to rewrite the field equations in terms of the dual potential u, 

{-d z u, d r Lo) = e n a{d r A, d z A). (25) 

As it was demonstrated in [8], using (25), we can rewrite (20), (21) in terms of a single complex 
differential equation with respect to the complex potential 

n 

£ = ea a + iuj, (26) 
In fact, we can go one step further and rewrite the field equations with respect to £. We get 

Re(£)Aa = -2Ae 2p a^v (27) 

■(<■*<) = S+*< (»> 

^■(cfi*^ = 0, n = 0...D-4 (29) 

2 ^-^ = H + ^ + ^ 4 *^ {u ~ v) (3o) 

(31) 

where we have redefined for convenience 2u = 2v — \ In a — ^. Equation (28) extends the Ernst 
equation [16] to the presence of a cosmological constant-indeed the extra term in Aa drops out 
if A = and we get (15). 

We stress here that Weyl coordinates (3) cannot be used once A / 0. We cannot therefore 
set a = R and the equations (27-30) are no longer integrable. We noticed however, that if we 
consider (16) then the Ernst potential of the Schwarzschild black hole (9) is simply £ = x in 
spheroidal coordinates. Given that we have the extension of Ernst's equation for A/Owe can 



now consider doing the same trick for A / 0. Consider first the four dimensional Kottler black 
hole with line element 



ds 2 = r 2i 



r 2 V(r) 



+ d9 2 ) - V{r)dt 2 + r 2 sin 2 9dij; 2 



dr 
rV(r) 



(32) 



dv 



and V(r) = 1 + k 2 r 2 — 2 ^-? Note that we will bypass Weyl coordinates by setting 
and not explicitly do the integral which involves elliptic functions. Comparing with (18) we 
now read off the metric components in (r, 9), 



a = r sin 



iW, £ = V, e 2 " = rV/ 2 . 



and therefore we set 



x- 1 



V{r) => x = -1 



X + l - k 2 r 2_2M 

r 

On the other hand setting as before y = cos 6 the metric (32) takes the form 



(33) 
(34) 



ds 2 



X 



1 



x + l 



dt l + r 1 



dx 2 



1 



+ 



dy 2 



(x 2 -l)(x + l) 2 ( r 2 k 2 + 2 ^-) 2 1-y 



+ (l-y 2 )# 2 



(35) 



and r is the real positive root of the third order polynomial, k 2 r 3 + — 2M = with 

respect to the new radial coordinate x. Note from (34) that there is an additional coordinate 
singularity when V = 1. We will come back to this for the more feasible case of D = 5. If we 
set k = we recover the usual spheroidal patch of (32) [3]. It is of more interest here to set 
M = in order to study adS and dS. For the former case we start with the global patch of adS 
and setting X = —x we get, 



ds z 



X + l 



dt 2 + 



dX 2 



+ 



X-l 2k 2 (X + 1){X - I) 2 X-l 



dy 2 



1-y 



+ (1 - y 2 )dcf 



(36) 



The coordinate range is X > 1, — 1 < y < 1 and only covers half of anti de-Sitter space with 
the boundary sitting at X = 1. The coordinate transformation for planar or hyperbolic slicings 
of adS goes through the same way. For A = we saw that in the (R, Z)-plane the integral 
curves x = constant are hyperbolas. So what happens for adS? To see this set kr = cosh ku 
which takes us to the usual global adS patch. Note now from the form of (1) that the Weyl 
coordinate r reads, 

r = ln|tanh-^| (37) 

which means that after a conformal transformation according to (8) we have X = coth(fcit) 
and y = cos(9). The integral curves are therefore similar but as we noticed earlier we need 
two patches to cover u G] — oo,+oo[ with a branch cut at u = 0. The adS boundary is then 



-Vdf 



+ 4h,+r 2 dK 2 D _ 2 



In the presence of a negative cosmological constant the generalised metric reads, ds^ — , , , , , 

with V(r) = k + k 2 r 2 — ^ a , adS curvature scale 2A = — (D — 1)(D — 2)k 2 and constant curvature of the D — 2 
compact space k = 0, 1, — 1. Note that in the presence of a negative cosmological constant we can obtain a black 
hole geometry with a flat (compact) horizon-the cosmological constant providing the necessary curvature scale. 
Furthermore, switching-off M gives us different slicings of constant curvature spacetime. 



at X = 1 and at X = — 1 for the second patch. For de-Sitter set k 2 = —a 2 starting from the 
locally static patch (32) to get: 



ds 2 



-dt 2 + 



dx 2 



+ 



x + l 2a 2 {x + l)(x 2 - 1) 1 + x \1 



dy 2 



+ (1 - y 2 )#< 



(38) 



The coordinate range is x > — 1 and the usual horizon at r = 1/a is not in this range. Therefore 
the spheroidals cover a yet smaller patch than the usual static de-Sitter (9). 
For D = 5 there is an interesting twist. Proceeding as before we get 



£ = V(r)(r cos 9) 



(39) 



and now £ is a function of r and 9. Note however that £q = r cos 9 is simply the Ernst potential 
for flat spacetime (i.e. when we set k = and M = 0) and we can therefore neglect it keeping 
the Ernst potential modulo flat spacetime (this also agrees with [18] for D = 5 and A = 0). In 
other words we once more have (34) (for D = 5) which yields, 



x + l 



£.2 r 2 _ 2M 



(40) 



Note that there is an extra coordinate singularity at r| = & (or 7 = 1) which is to the right 
of the event horizon at r = rjj or x = 1. This singularity goes to asymptotic infinity as the 
cosmological constant goes to zero k — > and dissapears for the flat case. Close inspection 
therefore reveals that for r > th there are two branches in spheroidal coordinates, one for 
ru < r < r s - , x > 1 and one for r > rg+, x < — 1. Setting y = cos 29 the first branch reads 



ds 1 



x-1 
x + l 



dt 2 + 



v / l + 2Mk 2 (x+ 1) 2 -1 



& 2 (x + 1) 



dx 2 



dy 



4(1 - y 2 ) 



_4(x 2 - 1)(1 + 2Mk 2 (x + l) 2 ) 
+ (1 - y)c^ 2 + (y + 1)# 2 



(41) 



and describes the region close to the black hole event horizon. In particular the limit M — > is 
ill-defined in this branch whereas we have a smooth k — > limit. For the second branch setting 
as before X = —x (ie. X > 1) and using the same angular coordinate, 



4±^ + 3 £ 



2Mk 2 {X - l) 2 + 1 
fc 2 (X- 1) 



dX 2 



4(X 2 -1)(1 + 2M£; 2 (X-1) 2 

dy 2 



■+ 



4(1 -y 2 



+ (1 - y)# 2 + (y + l)dV 2 



(42) 



This branch has the correct limit of M — > giving the equivalent of (36) in D = 5 and describes 
the region close to the adS boundary. An interesting open question, regarding the extension of 
spheroidals we have undertaken here, is if they can give the rotating generalisation of (9) as so 
happens for A = [16], [18]. 

As observed in [8] all equations (28-30) apart from (27) are independent of D whereas of 
course the metric (18) depends on the dimension. Therefore starting from a D = 4 stationary 
and axisymmetric solution with A = we can construct an infinite number of (D + n)- 
dimensional solutions parametrised by the n extra Weyl potentials ^-solutions of (29). In 
fact we have 

2(^(D),u - V(D+i),u)-^r = ~4*S. («<-»«). (43) 



where the subscript D refers to the spacetime dimension for v. Setting a = vm),u ~ v (D+i),- 
we rewrite the above equation in terms of R and Z in Weyl coordinates: 



4c7 z = R tf z 

8a, R = R(* 2 R - * 2 Z ) 



(44) 



In [8] examples were given where one started from a higher dimensional solution and found its 
seed D = 4 solution. Here we will do the opposite and uplift solutions making use of the Weyl 
formalism [12]. Consider as our seed metric in D = 4 that of Kerr, 



ds\ 



A , 2 n \ 2 sin i 

= ~ (at — asm Odip) H 7 - 

p z ' p z 



(adt - (r 2 + a 2 )dip)' 



+ f 



where M is the black hole mass, a the angular momentum parameter and 

A = r 2 + a 2 - 2Mr, p 2 = r 2 + a 2 cos 2 8. 
Comparing with (18) we now read off the relevant components 

^- p V/ 2 Va, 



O: 



A 



smtl, e 
2Ma sin 2 6 
p 2 -2M ' 



(p 2 - 2M) 2 
A/ sin 2 6 



(45) 

(46) 

(47) 
(48) 
(49) 



In order to uplift the solutions we will use Weyl coordinates (3) which are related to (45) by 
the relations 

R=\ r Ksm9 = a, Z = (r — M) cos 6. (50) 

Note that if we set A = ty/V6 in 29 and x = -4a/3 in (44) we obtain the 4 dimensional 
Weyl equations, (4) and (5) which are solutions of the static and axially symmetric metric (3). 
We have thus demonstrated the following: for each seed solution of A = 0, (27-30), for our 
example here (45), we can take any Weyl solution in 4 dimensions (in the form (3)), use the 
coordinate transform relating it to the seed coordinate system (here (50)) and thus obtain a 
D = 5 solution given by, 



ds 2 = p 2 K~^e^ 1 (dr 2 + dz 2 



+ R 2/3 e -A(R) 



-ez(dt + Adip) 2 + e~dip 2 



+ R 2/3 e 2A(R) # 2 

(51) 

We emphasize that R and Z are functions of r and z ie verify (50). 

Lets consider a flat potential (6) as an example. Start by noting that (4) has the obvious 
symmetry, A — > IX, I G R. As it was pointed out in [8], generically, when uplifting solutions 
asymptotic flatness is not guaranteed. Indeed in 5 dimensions, the extra dimensional dtp 2 
component in (51) reads R 2 / 3 e 2A and given the form of R in (44) an obvious problem will be 
to get this extra direction asymptotically flat. So lets start with A = HnR where I G R. It is 
straight forward to integrate and we get 



ds 2 = p 2 R^(dr 2 + dz 2 ) + R' 



-e~ (dt + Adtp) 2 + e 



2 dip 



+ R3 



2 + 2l di> 2 



(52) 



This is a solution to Einstein's equations for all I where one uses (50). For example, choosing 
/ = — 1/3 we get R 2 / 3 e 2A = 1. It is then easy to check that this is the Kerr string solution in 5 
dimensions. One can follow the same technique using this time the Weyl potential for the static 
black hole, (7) and adjusting the free constants in such a way as to obtain an asymptotically 
flat spacetime. It would be particularly interesting to extend this method to the case where 
A / 0. By analogy to what we have done here, it is possible that a starting point for such a 
generalisation would be to find how the additional fields of (27-30) transform upon keeping the 
Ernst potential £ fixed. 

Acknowledgments 

It is a great pleasure to thank my collaborators in [8] David Langlois, Daniele Steer and Dr. 
Robin Zegers as well as the organisers of the XBII conference. I also thank Roberto Emparan 
and Ruth Gregory for discussions on spheroidals. I finally thank the Galileo Galilei Institute 
for Theoretical Physics for their hospitality and the INFN for partial financial support. 

References 

[1] R. C. Myers and M. J. Perry, Annals Phys. 172, 304 (1986). G. T. Horowitz and A. Strominger, Nucl. 
Phys. B 360, 197 (1991). 

[2] R. Emparan and H. S. Reall, arXiv:hep-th/0608012. R. Emparan and H. S. Reall, Phys. Rev. Lett. 88, 
101101 (2002) [arXiv:hep-th/0110260]. H. Elvang, R. Emparan, D. Mateos and H. S. Reall, Phys. Rev. 
Lett. 93, 211302 (2004) [arXiv:hep-th/0407065]. 

[3] T. Harmark, hep-th/0408141. Phys. Rev. D 70, 124002 (2004) [arXiv:hep-th/0408141]. 

[4] R. Emparan and H.S. Reall, Phys. Rev. D65 (2002) 084025, hep-th/0110258. 

[5] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] [arXiv:hep- 
th/9711200]. 

[6] S.W. Hawking, C.J. Hunter and M.M. Taylor-Robinson, , hep-th/9811056 G. W. Gibbons, H. Lu, D. N. Page 
and C. N. Pope, Phys. Rev. Lett. 93, 171102 (2004) [arXiv:hep-th/0409155]. G. W. Gibbons, H. Lu, 
D. N. Page and C. N. Pope, J. Geom. Phys. 53, 49 (2005) [arXiv:hep-th/0404008]. 

[7] C. Charmousis and R. Gregory, Class. Quant. Grav. 21, 527 (2004) [arXiv:gr-qc/0306069]. 

[8] C. Charmousis, D. Langlois, D. Steer and R. Zegers, arXiv:gr-qc/0610091. 

[9] T. Tanaka, Prog. Theor. Phys. Suppl. 148, 307 (2003) [arXiv:gr-qc/0203082]. R. Emparan, A. Fabbri and 
N. Kaloper, JHEP 0208, 043 (2002) [arXiv:hep-th/0206155]. 
[10] A. L. Fitzpatrick, L. Randall and T. Wiseman, arXiv:hep-th/0608208. 

[11] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers and E. Herlt, , Cambridge University Press, 1984. 
[12] H. Weyl, Annalen Physik 54, 117 (1917). 
[13] H. Quevedo, Phys. Rev. D 33, 324 (1986). 

[14] W. Israel and K. A. Khan, Nuovo Cim. 33, 331 (1964). R. C. Myers, Phys. Rev. D 35, 455 (1987). 

H. F. Dowker and S. N. Thambyahpillai, Class. Quant. Grav. 20, 127 (2003) [arXiv:gr-qc/0105044]. 
[15] T.Lcwis, Proc Roy Soc Lond. A 136, 176 (1932). A Papapetrou, Ann. Phys. 12 (1953) 309. A. Papapetrou, 

Ann. Inst. H. Poincare A4 (1966) 83. 
[16] F.J. Ernst, Phys. Rev. 167 (1968) 1175. F.J. Ernst, Phys. Rev. 168 (1968) 1415 

[17] P. T. Chrusciel, G. M. Greuel, R. Meinel and S. J. Szybka, Class. Quant. Grav. 23, 4399 (2006) [arXiv:gr- 
qc/0603041]. G. Neugebauer and R. Meinel, "Progress in relativistic gravitational theory using the inverse 
scattering J. Math. Phys. 44, 3407 (2003) [arXiv:gr-qc/0304086]. R. Meinel and G. Neugebauer, Class. 
Quant. Grav. 12, 2045 (1995) [arXiv:gr-qc/0302114]. 

[18] T. Mishima and H. Iguchi, hep-th/0504018 

[19] D.M. Zipoy, J. Maths. Phys 7 (1966) 1137. 

[20] F. Kottler, Ann. Phys. 56 401. 



